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A sharp maximum in the critical current J c as a function of temperature just below the melting point of the Abrikosov flux 
lattice has recently been observed in both low and high temperature superconductors. This peak effect is strongest in twinned 
crystals for fields aligned with the twin planes. We propose that this peak signals the breakdown of the collective pinning 
regime and the crossover to strong pinning of single vortices on the twin boundaries. This crossover is very sharp and can 
account for the steep drop of the differential resistivity observed in experiments. 
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The discovery of the high-temperature copper-oxide 
superconductors has renewed the experimental and the- 
oretical interest in the properties of the mixed state of 
type-II superconductors in a magnetic field Vari- 
ous experimental techniques, including standard current 
versus voltage curves, are used to measure the critical 
current density J c needed to depin the flux-line array 
and to investigate its temperature and field dependence. 
Naively J c is expected to decrease monotonically as the 
temperature or the applied field are raised towards the 
mean field H C 2(T) line. It has, however, been known 
for some time that an abrupt increase in J c as a func- 
tion of field or temperature can occur in conventional 
low temperature superconductors near H C 2 [||. A qual- 
itative explanation of this phenomenon, referred to as 
"peak effect" , was proposed a long time ago by Pippard 
H , who argued that the increase in J c is associated with 
the softening of the shear modulus cqq. A more quanti- 
tative explanation of the peak effect as arising from the 
softening of all the elastic moduli of the flux lattice near 
H C 2 was presented by Larkin and Ovchinnikov Q . 

Recently a sharp maximum in J c as a function of tem- 
perature has been observed in both untwinned ]^-|^] and 
twinned || YBCO crystals, as well as in some low tem- 
perature superconductors 



10 11 1 . The new feature is 
that in this case the peak occurs away from just 
below the temperature T m where the flux lattice melts 
into a flux-line liquid. In view of the old suggestion by 
Pippard || and the work by Larkin and Ovchinnikov j| , 
it is natural to associate it with the softening of the shear 
modulus C66 that occurs at the melting point. In twinned 
YBCO crystals the peak depends strongly on the orien- 
tation of the applied field relative to the twin planes: it is 
largest for flux motion along the twin planes and external 
fields aligned with the c axis and it weakens as the field 
is tilted out of the plane of the twins In untwinned 
YBCO single crystals the peak is much smaller: it shifts 
towards T c and becomes less pronounced as the sample 



purity is increased ||. These observations suggest that 
strong anisotropic pinning centers such as twin bound- 
aries enhance the peak effect. 

In this paper we first examine in detail the temper- 
ature dependence of the critical current from collective 
pinning by point defects |f^,[L3|,|lj , and show that J c can 
exhibit a sharp rise near T m for a narrow range of mag- 
netic fields due to the abrupt decrease of the shear mod- 
ulus. This may provide a mechanism for the small peak 
effect observed in untwinned single crystals. New results 
on anisotropic collective pinning in samples with a family 
of parallel twin planes are also presented and show that 
the same mechanism is in principle operative in twinned 
samples, when vortices are pinned collectively by an ar- 
ray of twin planes. On the other hand, collective pinning 
is very weak in this case and cannot account for the large 
increase in J c observed in these samples. In the second 
part of the paper we show that in twinned samples a large 
peak in J c near melting can arise from strong pinning of 
individual vortices on the twins. This mechanism is en- 
hanced at T m by the vanishing of cqq and can account for 
the sharp drop in the resistivity observed in experiments. 

The critical current density of an elastic medium 
pinned by weak disorder can be calculated using the col- 
lective pinning theory jjj. Weak disorder destroys the 
translational order of the flux-line lattice and results in 
the coherent pinning of vortex bundles of extent L c and 
R c in the directions parallel and perpendicular to the ap- 
plied field, H. The pinning lengths R c and L c , defined 
as the distances at which the lattice distortion due to 
disorder is of the order of the range £ of the pinning po- 
tential, are determined in terms of the elastic constants 
of the lattice by balancing the elastic deformation en- 
ergy against the pinning energy. The critical current J c 
is then defined as the current where the Lorcntz force 
balances the pinning force, or BJ c /c « y/W/V c , where 
W — n p < f 2 > is the mean square pinning force, with 
n p the volume density of pins and / the elementary pin- 
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ning force, and V c — R 2 L C . 

We consider a three dimensional flux-line array in a 
sample with an external magnetic field aligned with the 
c axis, which is chosen as the z direction. Disorder 
is described as a quenched random potential per unit 
length V(r) with zero mean and Gaussian correlations, 
V(r)V(r') — r(r, r'). The overbar denotes the disor- 
der average and the correlator T(r, r') is determined by 
the strength and geometry of the disorder. The static 
elastic deformation of the lattice due to disorder can 
be evaluated by a perturbation theory in the pinning 
potential |12,Bj. To lowest order in perturbation the- 
ory the components of the mean square displacement 
Uij(r) =< Aui(r)Auj(r) > induced by the random po- 
tential (here Attj(r) = Uj(r) — itj(O) and the brackets 



denote a thermal average) are given by jlj , 
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(2tt) 3 J (2tt) 3 
G lk (q,LU = 0)G jl (rf,W = 0)W kl (q + q'), (1) 



with Gij(r,t) the elastic Green function of the lattice, 
f p (r) = — n(r, t) VV r (r) the pinning force per unit vol- 
ume and Wfc/(q,q') = /p,i(q)/pj(q') the pinning force 
correlator. Here n(r,t) = J2 n ^ 2 H r i- ~ r n(z,t)), with 
r = (rj_ , z) , is the coarse-grained microscopic vortex den- 
sity field (the flux lines are parametrized by their trajec- 
tories {r n (z,i}}). The main contribution to Eq. ([[]) for 
the case of interest below comes from the transverse part 
of the elastic Green's function, given by 



Gy(q,w) 



-iuj( + c 66 qj_ 2 + c 4 4(qj_, q z )q 2 
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where = <% - q±iq±j, with q ±i = qu/qx, and c 66 
and C44(<7j_, q z ) are the shear and tilt moduli of the vor- 
tex lattice, respectively. Thermal fluctuations can be in- 
corporated in the perturbation theory by separating out 
in the vortex positions r n (z,t) the deviation from equi- 
librium due to pinning from that due to thermal effects, 
as described in The main effect of thermal fluctu- 

ation is the replacement of the upper cutoff <?o ~ in 
the wavevector integral by a thermal cutoff qr ~ (£ 2 + < 
u 2 > th y 1/2 = qoil+T/Tdp)- 1 / 2 . The depinning tem- 
perature T dp is defined by < u 2 (T dp ) > t h~ £ 2 , where 



W lj {q) = WS lJ (2TT) 3 S^{q), with W = l/^a 2 . The 
pinning lengths and the critical current for this disor- 
der geometry have been calculated elsewhere Q, but 
it is instructive to display the dependence of J c on 
the elastic constants. For low defect densities {R c > 
A - this is the so-called large bundles regime), the 
dispersion of the tilt modulus can be neglected and 
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Here 



jo = cH c /3^nX is the depairing current and j sv — 

JoiWa^/e^) 2 3 is the single vortex critical current den- 
sity. The critical current contains both an explicit tem- 
perature dependence from the thermal smoothing of the 
pinning potential and an implicit T dependence through 
the superconductor's parameters that determine the elas- 
tic constants. To extract the strong temperature depen- 
dence of C66 near melting, we write Cqq = Cg 6 r(T/T m ), 



where 



66 



depends weakly on T and the function 



drops sharply from unity to zero at the clean flux lat- 
tice melting temperature T m . This is defined by < 



u 2 (T m ) >th~ c 
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with cl ~ 0.1 — 0.3 the Lindcmann 



parameter. The critical current is then J c ~ -j 
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temperature dependence of J c near T m is controlled by 
the parameter a = T m /T^ p = (c L a /£) 2 . For H = 6T 
and temperatures near T c , we use £ w f 00 A and cl ~ 0.2 
to obtain T m /T^ p « 0.1. At low temperatures (T « 
T dp ,T m ) Jc decreases very slowly with T. At higher 
temperatures, but still well below T m , the elastic con- 
stants are only very weakly temperature dependent and 
the temperature dependence of J c is controlled by ther- 
mal fluctuations, yielding a decrease of J c as T grows. 
As T m is approached from below, cqq softens and the 
flux lattice can better adjust to the pinning centers, rais- 
ing J c . Finally, near T m the function r(T) and there- 
fore T dp drop sharply in a narrow temperature range giv- 
ing rise to an increase of the pinning barriers and an 
associated drop in J c . For larger defect densities the 
dispersion of C44 is important (this is referred to as the 

small-bundle regime). In this case J c jo 6 ^^ C66 ( 1 + 



T/T dp ^j exp 
Again, 
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the critical current can be written 



T/T dp ) 



< u 2 (T) > th is the mean square thermal excursion of exp{— ar 3 / 2 [l + T/(T$ p r 1 / 2 )] 3 }, where a is practically in 
the vortices about their equilibrium positions. When the 
vortex array is described as an elastic continuum, T dp = 



2a £ 2 v / c 6 6C44 = 2ao£ 2 v / C66C44/A, with a = \/4>o/B the 
mean intervortex separation, A the penetration length in 
the ab plane and C44 = ca{q\_ = 0, q z = 0) Q. 

When the flux array is pinned by isotropic point dis- 
order, the random potential is short-ranged in all direc- 
tions, with T(r, r') = 7<S( 3 )(r-r') and 7 w (U tl 3 ) 2 n p [l + 
C(n p £ 3 )] , with Uq the depth of an individual pinning 
potential. The pinning force correlator is isotropic, 



dp' 

dependent of temperature near T m . We find J c ~ e~ a 
below T m and J c ~ e" a a as T -> T~, yielding a sharp 
rise of J c in a very narrow temperature range. This mech- 
anism may be responsible for the small peak effect ob- 
served in untwinned single crystals. We remark, however, 
that the melting transition is broadened by the presence 
of defects and the drop in critical current is more likely 
associated with the onset of plastic motion of vortices. 

We now consider collective pinning in a sample with 
a single family of twin boundaries of mean separation d 
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spanning the zy plane for fields aligned with the twin 
planes. This is the experimental geometry where the 
peak in J c is strongest |1. We define two different pin- 
ning lengths, i? c j| and R c ±, corresponding to the size of 
the vortex bundle in the directions parallel and transverse 
to the twin planes, respectively. If R c ± >> d, pinning oc- 
curs via the collective action of many twin planes. Each 
twin is described as a sheet with a large concentration of 
point defects. The correlator of the random potential is 
given by T(r, r') = "i\g{\x — x'\)S(y — y')S(z — z'), where 
7i ~ {Uo^ 3 ) 2 n^ [l + 0(£/d)] is proportional to the areal 

(2) 

density n p of pins on each twin plane, and g(x) describes 
correlations in the distribution of twin planes. On dis- 
tances large compared to the twin spacing d the twins are 
essentially uncorrelated [g(x) w (1/ 'd)d(x)] and the pin- 
ning force correlator is Wij(q+q') = WrSij (27r) 3 5^ (q+ 
q'), with Wt = 7i/£ « a o- Collective pinning in this 
regime is very similar to collective pinning by point de- 
fects in bulk. The dependence of J c on the elastic con- 
stants and temperature is identical to that obtained for 
isotropic point disorder, with the replacement W — ► Wt- 
A peak effect in densely twinned samples may then in 
principle arise from the same mechanism discussed above 
for untwinned crystals. On the other hand, the mean 
squared pinning force Wt is still determined by the ef- 
fective volume density of pins, which is now given by 
rip ao/^d. As a result, the anisotropy due to the twin 
planes increases the pinning volume of a factor (d/d p ) 3 , 
with dp ~ (rip)" 1 / 3 , correspondingly decreasing the criti- 
cal current. For this reason collective pinning in twinned 
samples is very weak, especially if d >> ao and cannot 
account for the observed critical currrents. 

The dominant pinning mechanism in twinned crystals, 
particularly in sparsely twinned samples, is the strong 
pinning of individual vortex lines trapped on the twin 
boundaries. As T m is approached fom below intervor- 
tex interactions weaken and the vortices on the twins 
become more strongly pinned than those in the chan- 
nels between twins. The main contribution to the criti- 
cal current arises then from pinning of single vortices on 
the twins, and particularly from those vortex segments 
that are strongly pinned in rare regions with an excess 
of impurities. As described below, it is the rise in the 
fraction of such strongly pinned vortex segments on the 
twins with T that can be responsible for the peak in J c 
in twinned samples. 

To evaluate the critical current due to strong pinning 
in regions with excess impurities, we consider a repre- 
sentative vortex line trapped near a twin plane and in- 
teracting with its neighbors at an average distance ao in 
the lattice. The remainder of the lattice, even though 
not directly pinned by the twin, is held in place by in- 
teractions. The magnitude of the elastic force associated 
with displacing a length L of the representative fluxon a 
transverse distance u from its equilibrium position is 



F e i(u,L) ~ ?i— +C6 6 uL. (3) 

The first term is the force associated with tilting the 
representative vortex, with l\ the tilt energy per unit 
length. The second term arises from the interaction with 
the neighbors. The typical pinning force exerted on a 
vortex segment of length L is (FjIT)) 1 / 2 - (Wx^L) 1 / 2 , 
with Wi — 7i/£ 5 ao the mean squared pinning force per 
unit volume due to a single twin of thickness £. The most 
effective pinning arises from rare regions with an anoma- 
lously large impurity concentration that pin strongly the 
vortex segment. The pinning forces F p in these regions 
exceed the typical pinning force {F 2 (L)f/ 2 and give the 
dominant contribution to J c . The problem of strong pin- 
ning of vortex lines is analogue to that of incommensurate 
charge density waves and can be rigorously discussed fol- 
lowing Ref. | jl7| |. Here we prefer, however, to follow the 
more phenomenological, but physically intuitive discus- 
sion given by Coppersmith |Tq| . The condition for the 
strong pinning is F p (£) > F e t(£) fl^9|. The critical cur- 
rent J c is proportional to the density n of the strongly 
pinned vortex segments where this condition is satisfied. 
To find it we note that the pinning force F p scales as 
the impurity excess in the region and it can be shown to 
be Gaussian-distributed with variance {F 2 ) 1 / 2 |^8| . The 
density n of vortex segments that are strongly pinned in 
these excess-impurity regions is then, 

poo poo 

n~ / dL dF p e-^/^ F i^\ (4) 

Jo JFel(L) 

where we have used u ~ £. Using Eq. [|, one can show 
that the integral over L is dominated by the length scale 
L* ~ y^x] c 66 where the single-vortex tilting force bal- 
ances the force of interaction of the pinned vortex with 
the rest of the lattice. The critical current J c is then 
given by, 

J c ~ n ~ exp[-(c 66 /c^) 3/2 ], (5) 

where c* m = [2W\j \fT\) 2 ^ . As T approaches T m , the 
shear modulus softens and drops to zero. Correspond- 
ingly, J c grows according to Eq. (^). The result given 
in Eq. (||) applies for cqq > Cg 6 . The mechanism of 
strong pinning just described is also operative in un- 
twinned samples where vortices are pinned by isotropic 
point disorder. We note that strong pinning of isolated 
vortices yields the same functional dependence of J c on 
C66 as collective pinning of small bundles. In both cases 
J c ~ exp(— aegg 2 ). In untwinned crystals a is essen- 
tially the same for these two pinning mechanisms. In 
twinned crystals collective pinning of vortex bundles is 
very weak (a ~ l/W') and strong single vortex pin- 
ning in regions with excess impurities (a ~ 1/Wx) con- 
trols the critical current even when the condition for 
single-vortex pinning is satisfied only locally. In this 
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case the main contribution to J c arises from pinning 
energy barriers which are large compared to the typi- 
cal barrier E p (L*) — (W / iC 2 /c5 6 )(cg 6 /c 6 6) 1 ^ 4 , but small 
compared to the scale of the elastic energy of interac- 
tion of the vortex segment with the rest of the lattice 

Most experiments do not measure the critical current 
as defined theoretically, but rather the nonlinear resis- 
tivity. For comparison with experiments it is important 
to discuss the small finite resistivity due to the creep of 
strongly pinned vortex segments at low measuring cur- 
rents, J < J c - To find the energy barriers that determine 
the creep rate, we consider the distribution of barriers 
separating different metastable states in the regime of 
strong pinning. As discussed above for pinning forces, it 
has been shown that the impurity energies 5E of strongly 
pinned vortex segments of length L are described by 
a Gaussian distribution, - exp[-(<S£') 2 /2H/ r i J L£ 4 ] @. 
A gain SE in pinning energy is associated with a cost 
8T e i = ii£_ 2 /L + cqq^L in elastic energy. The result- 
ing activation barrier is therefore E p = SE — 5J- e i . The 
distribution function of creep barriers can be evaluated 
following H], 

P(E p ) ~ / dLe-^+s^) 2 /2W lL e „ e -^7^. ( 6 ) 
Jo 

The typical creep barrier at J ~ J c is E p ~ W\€,/cqq, 
At arbitrary currents the creep barrier becomes E P (J) ~ 
Wi^/cQef(J/J c ). The function f(J/J c ) decreases as J/J c 
grows, but its explicit form cannot be obtained by our 
dimensional analysis. The resulting thermally activated 
resistivity is 

with J c given by Eq. (|^). The temperature dependence 
of the resistivity at a fixed current is governed by the 
temperature dependence of J c . When cqq decreases near 
melting, the creep resistivity drops exponentially. This 
result applies for J < J c and temperatures near, but 
below melting, where Cqq is not too small, Cqq > \/Wi£. 

In the regime just described dominated by strong pin- 
ning of individual vortex segments on twin planes, in- 
teractions are still strong enough to hold the lattice to- 
gether so that the remainder of the flux array can be 
described as an elastic continuum. On the other hand, 
in the presence of an applied current the competition 
between strong single- vortex pinning at the twin bound- 
aries and the elastic deformations of the portions of flux 
lattice between twin planes eventually leads to the devel- 
opment of large strains in the regions next to the twins. 
This results in the break-up of the lattice and the on- 
set of plastic flow. The condition for the onset of plastic 
flow is obtained by considering the elastic deformation of 
a flux lattice in the channel < x < d between two twin 



planes driven by a current transverse to the twins, yield- 
ing a Lorentz force fi, — BJ/c per unit volume in the 
+y direction. The displacement u y (x) of the flux lattice 
is the solution of C,d t u y — c^d^Uy + f^, with boundary 
conditions u y (0) — u y (L) = 0. In an elastic continuum 
the mean velocity v y = dtu y must be spatially homoge- 
neous and u y (x) = (/i/2c66)^(d — x). The lattice yields 
when the strains at the twin boundaries become too large 
and the condition \d x u y (x)\ « 1 is no longer satisfied. 
An approximate condition for the onset of plastic flow 
can then be written as f^d ~ cqq, with a corresponding 
current scale, J p ~ c§§jBd. The onset of plastic flow 
at J ~ J p corresponds to a sharp rise in the differential 
resistivity and the critical current. 
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